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Abstract 

Rigorous UTD (Uniform Geometrical Theory of Diffraction) diffrac- 
tion coefficients are presented for a coated convex cylinder simulated 
with generalized impedance boundary conditions. In particular, ray 
solutions are obtained which remain valid in the transition region and 
reduce uniformly to those in the deep lit and shadow regions. These 
involve new transition functions in place of the usual Fock-type inte- 
grals, characteristic to the impedance cylinder. A uniform asymptotic 
solution is also presented for observations in the close vicinity of the 
cylinder. As usual, the diffraction coefficients for the convex cylinder 
are obtained via a generalization of the coirepsonding ones for the 
circular cylinder. 
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I Introduction 

The problem of scattering by a smooth convex impedance cylinder has re- 
ceived much attention. Wang [1, 2] presented ray-optical solutions for the 
impedance and coated cylinders. His results are valid only in the deep 
lit and shadow regions and do not apply to the case where the observation 
point is in the transition region. Wait and Conda [3, 4] developed a solution 
which is valid in the transition region and for observation points on and 
off the surface. However, as pointed out by Pathak [5] it did not uniformly 
reduce to the ray solution [6, 7] exterior to the transition regions. Also, it 
is not valid on the portion of the surface in the transition region and these 
limitations were the primary motivation in Pathak’s work [5] for the per- 
fectly conducting convex cylinder. Recently, Kim and Wang [8] presented 
a solution applicable to a coated cylinder that remained valid in the tran- 
sition region. They employed a heuristic approach to obtain the numerical 
values of the resulting transition integral applicable to a coated cylinder. 
Their solution is uniform but is not applicable to the close vicinity of the 
cylinder. 

In this paper we present a rigorous UTD solution of the diffraction by 
a coated cylinder simulated with generalized impedance boundary condi- 
tions. In addition, a uniform asymptotic solution is obtained which remains 
valid when the observation point is in the close vicinity of the cylinder. An 
important aspect of the paper is also the use of second order generalized 
impedance boundary conditions (GIBC) for the simulation of the coating. 
Their derivation has already been given in [9] and [10] and are character- 
ized by the inclusion of higher order field derivatives in their definition. 
Because of this they are less local which leads to an improved simulation 
(with respect to the standard impedance boundary condition - SIBC) of the 
coating in a manner analogous to the order of the highest derivative kept 
in the condition. Recently, they were successfully applied to a number of 
diffraction problems [11], [12] and have also been used in numerical simula- 
tions of multilayer coatings (see fig. 1) [13]. These applications provided a 
measure of the accuracy of the proposed GIBC and in particular accuracy 
criteria were derived in [13] for the second order conditions as a function 
of coating thickness and composition. 

The UTD solution to be presented here parallels that given by Pathak [5] 
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for the circular perfectly conducting cylinder. However, in the case of the 
coated cylinder the resulting UTD expressions are in terms of Fock-type 
integrals whose efficient evaluation is of primary interest. In the follow- 
ing we first present the eigenfunction solution based on the second order 
GIBC simulation of a circular coated cylinder. By employing Watson’s 
transformation this is written in integral form which is then cast in a ray 
representation. The ray solution is subsequently generalized to the case of a 
general convex cylinder. Finally, the evaluation of the Fock-type integrals 
is discussed and some results are presented which validate the accuracy 
of the GIBC eigenfunction and ray solutions. In the process, we demon- 
strate the improved accuracy of the GIBC solution over the corresponding 
SIBC solution, and it is also shown how the presented UTD solution can 
be extended to treat multilayered coated cylinders. 


II Boundary Conditions 


Consider the coated circular cylinder shown in figure 2(a). We propose to 
simulate the cylinder with an equivalent one (see fig. 2(b)) satisfying the 
boundary conditions (an e JU,< time convention is assumed throughout) [13] 


_ (Q2 + Qq) 7 TT &Hz 

* a x ZHz k*Va x dp 

for H z -incidence or 

E _ (<*2 + Qq) „ jj q 2 &Ex 

z a, * kW(a' 2 + a' 0 ) dp 


( 1 ) 

( 2 ) 


in the case of ^-incidence. In these Z = 1/Y is the free space intrinsic 
impedance, k = y- is the free space wave number, and a m with a' m are 
constants which are dependent on the material properties of the coating. 

Assuming a single layer coating of thickness 6 and denoting the associ- 
ated relative constitutive parameters as e T and /x r , we have 

a c = kS (3a) 


q i = ~j 


(36) 


3 



a 2 = kS/tr 


and 


(3c) 

(4a) 
(46) 
(4c) 

These constants are valid for ^/c r /i r 6 ~ 0.15A and for higher contrast coat- 
ing they must be redefined as [9] 


a'=l 


a\ = +jkfi r 6 


a 2 ~ 0 


= { N - 577) taE(HJV) - tan (|£)] 

a x = —je T |l + tan(k6N ) tan 

i[M«»0-t-(^)+H(w-i) 

1 1 + tan (k6N) tan j Jj 


a 2 = 


a' = 

0 \ 2 N) 


1 + tan (kSN) tan 


(k6_ 

\2N 


)] 


a\ = jfi tan (kSN) - tan 


(In 


2N 


1 + 


taa(kSN) tan (jf) ~ w ( N " 


■ |tan(WJV) - tan f ) 


(5a) 

(56) 


(5c) 


(ба) 

( бб ) 


(6c) 
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We note that (1) and (2) reduce to the standard impedance boundary 
conditions (SIBC) when 02 and a ' 2 are set to zero and 


o >0 a x 
V ~ a’ 


( 7 ) 


where r) = —j tan( £ yJt T [i T S ) is the normalized impedance of the coating. 

We also remark that a m and a' m satisfy the duality relation 


q i = Q 2 + < 

a 2 + a 0 a\ K ’ 

The aforementioned boundary conditions can be employed for simulat- 
ing multilayer coating by simply redefining the material constants a m and 
a' m as discussed, for example, in [10] and [13] 


III Eigenfunction and Integral representa- 
tion 

Consider the plane wave 

tx’ = u 0 e ikx = u 0 e’ kpco *+ (9) 

incident on the coated circular cylinder shown in figure 2, where u 0 is the 
amplitude of the wave and u' denotes either the E z or H z component of 
the incident field. An eigenfunction representation of the total field in the 
presence of the cylinder can be generally written as 

OO 

«'= £ r[uk(,)+A^\kp)\t-‘"* ( 10 ) 

n=— 00 

where J n (‘) is the Bessel function of order n and denotes the nth 

order Hankel function of the second kind. To find the constants A n we 
enforce the GIBC given by (1) or (2) at p — b. This yields 

J'Jkb) + Q(n)J n (kb) 

Hn^'(kb) + Q(n)Hn\kb) { ’ 
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in which the prime denotes differentiation with respect to the Bessel or 
Hankel function argument, 


f -3 “ n2 pf^r) for H*- incidence 


Q(n) = 


(12) 


[ —j - n 2 j„< ) for ^-incidence 

and when a 2 = a' 2 = 0, Q(n ) reduces to 

( —j 2 * = —jrj for ^-incidence 


Q(n) = Q = 


— j 2 ? = — j- for ^-incidence 


(13) 


To obtain a ray representation of (10) in conjunction with (12) we must 
first recast tt* in integral form. Employing the usual Watson’s transforma- 
tion [14] and noting that Q(n) = Q(— n), we find [5] 


u 1 = tii + u 2 


where 


J'„(kb) + Q(v)J v (kb) 




or 
u l 
and 


Uo !-J[ J ^ kp) H^'(kb) + Q(v)H (2 \kb) 

= * r k'(^) - 

2 •/-<»[ (B) + Q(v)Hl \kb) 


e-^di<14a) 


(146) 


l Hl 2) (kb) + Q(v)Hl 2 \kb) 


e ,»«l2 


e -jv(2r+4>) _|_ e -jv(2*-4 

1 - e-> 2,/ * 


(15) 


- dv 


in which 


7 r 


V> = \<t>\~ 2 


( 16 ) 
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and Q(v ) is given by (12) upon letting n — ► v. The field component tii, 
defined above, includes the geometrical optics and dominant surface diffrac- 
tion contributions whereas t *2 denotes the creeping waves which circulate 
around the cylinder more than once. We are interested in an asymptotic 
evaluation of u\ and to do this we must separately consider each of the 
space regions illustrated in fig. 3. 


IV Field in the Deep Lit and Shadow Re- 
gions 


In the deep lit region, the geometrical optics incident and reflected rays 
represent an accurate first-order high-frequency approximation of the total 
field. The geometrical optics field can be extracted from ui upon performing 
an asymptotic evaluation of (14) yielding [18] 


u ao (P L ) a u‘(P L ) + u'(Qr)R 


in which 


R, 

Rh 


lA \h + e 


- — c~ jW 


a' — a'j cos O' + a' 2 cos 2 6 l 
a' + a\ cos O' + a' 2 cos 2 6 l 

a 0 — ai cos O' + 02 cos 2 6' 
a„ + ai cos O' + 02 cos 2 


(E x - incidence) 
(H x - incidence) 


(17) 


(18a) 

(186) 


axe the reflection coefficients associated with the coating [9]. The parameter 
O' is illustrated in fig. 4 and t is the distance (must be large) from the surface 
reflection point Qr to the observation point Pi. Finally, p is the reflected 
ray caustic distance and for a convex cylinder it is given by 


- PgiQfi) cos & 
P = 2 


(19) 


P 9 (Qr) is the radius of curvature of the surface at Qr and is equal to b for 
the circular cylinder. 
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For observations in the deep shadow region, a residue series representa 
tion of is more appropriate. From (14), we obtain 


4 00 

“> = -“- 11 £ 


( 20 ) 


b! 2 (kb)& \H^\kb) + q(u)hP(U)] _ 

where i/ m are the zeros of the transcendental equation 

+ Q(v„)H<£m = 0 ( 21 ) 

Obviously, (20) does not provide a ray-picture interpretation of the creeping- 
wave diffraction which is desirable for generalizations to non-circular cylin- 
ders. To recast (20) in a form compatible with the Keller type GTD for- 
mat, the Debye approximation must be employed for the Hankel function 
for kp » \u m \. This yields 


ui(P.) * u\Q x )T tX 


D-jk* 


( 22 ) 


where 


T.J, = -£ ©„(<?!) • t-*"' • V m (Q,) 


(23) 


m=l 


is the corresponding diffraction coefficient and for the circular cylinder 
V m {Q x ) = V m (Q 2 ) 


(24) 


/Li 

Vtt kH 


e 1 * 


1 h 


[V * k^b H l 2 J(kb)l {Hl 2) (kb) + Q{u)H (2 \kb)} v=vm j 

with 0 as defined in figure 5. 

A generalization of (22) - (24) for convex cylinders is now possible by 
making the replacements 

b - p a (Q h2 ) (25) 


e~* Vm6 


exp 



Moo PiV ) . 


(26) 


Consequently, the attachment coefficient 'D Jn {Q \ ) is no longer equal to the 
launching coefficient T> m (Q 2 ). 
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V Field in the Transition Region 


The geometrical optics and creeping wave solutions presented above are not 
valid in the transition region close to the shadow boundaries as illustrated 
in figure 3. New uniform expressions are, therefore, required to overcome 
this limitation. These can be derived following the procedure adopted by 
Pathak [5] for a perfectly conducting convex cylinder. For the lit region we 
find 


«i (Pl) w u\Pl) + u t (Qn)R tt h\j-z 

where R St h is now given by 

Ra,h = 


+ e 



F(x) = 2jy/xe ’ x j e *** dx 


x' = —2 m(Qp) cos 8' 


™W*>- 


and 


P = 


P 3 (Qr) COS 8' 


r( > n\- e — T V ( T )-n T )nT) B - jX T, 

v/tF J- oo WI(t) - q(T)W 2 (T ) 


in which 


q( T ) = mQ( u) 


jkt 

(27) 

, e‘)} + G(x', q ) 

(28) 


(29) 


(30) 


(31) 


(32) 

is defined by 


■ jxr dT 

(33) 


(34) 
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( 35 ) 


whereas V(r) and W li2 (r) are the Fock-type Airy functions [15] 


2 jV(r) = Wi(t) - W 2 {t) 

W,. 2 (r) = -Lf e"- 1 ’' 3 * (36) 

yTT ^ri,a 


The contour Ti runs from ooe ^ to oo—je and r 2 is the complex conjugate 
of r 2 . We remark that in the case of an SIBC simulation, q(r) becomes a 

constant. 




For the shadow region, we have 




Ut(P t ) w 

V* 

(37) 

where the diffraction coefficient 

n: i 

is now redefined as 
g-if 


T s ,h = 

— \Jm(Qi)m(Q2)\j—e~i kt 

2x v^r l 1 ~ F(ksa)} + G(x, q ) 

(38) 

in which 





x = 

MQi) Pg{t') 

(39) 


t = / dt' 

Jt'(Qi) 

(40) 

and 


X 2 

(41) 


2m(Q 1 )m(Q 2 ) 


As is usually the case, (27)-(40) were first derived for the circular cylin- 
der and were subsequently generalized for non-circular convex cylinders. 

VI Field in the Close Vicinity of the Conex 
Cylinder 

In all of the above derivations we have assumed that l is large. Conse- 
quently, the given expressions are not adequate for field computations in 
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the near vicinity of the cylinder. In this case, it is possible to obtain a 
suitable approximation of the integral (14) by replacing the Hankel and 
Bessel functions in terms of Airy integrals and then employing a Taylor 
series expansions for these integrals. 

Following a procedure similar to the one given in [5], the resulting ex- 
pressions for an arbitrary smooth convex cylinder are 


u,(P) « u\P N ) 


r ihx ' - t, (—^(jhz'y + e=^ {A ,(z') - A 2 (z')} 


n= 0 


n! 


(42) 


when Pn is in the lit region, and 


u 1 (P)Ku i (Q 1 )e~ jkt 


£s(£W)l 6 

.Pg(Q i). 


(Ai(z) - A 2 (z)} 


(43) 


when Pn is in the shadow region. In (42) and (43) 

A ,(D) = 5,(D) + i^g’AD) - ~g,(D) - (D) - j*^g[(D) (44) 


A 2 (D) = hg,(D) + j | S ;(D) - 2^g,(D) - ^ S J(X>) (45) 


1 r°° 

= 7?/- 

1 f°° 


e ~}Dr 

W&T) - dTjm ~) dT 


q{T)e-*» 
W&t) - 


(46) 


(47) 


z' = —m(pN) cos 8' 
for Pn in the lit region and 


= 

MQi) Pai’ 1 ') 


(48) 


( 49 ) 
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for Pw in the shadow region 


t 


-L 


*'(Pn) 


dt' 


<'(<?!) 


kdjPs) 

m(P N ) 

d{Ps) = P ~ P 9 {Pn) 


(50) 

(51) 

(52) 


The points P and P ^ in the lit and shadow regions are illustrated in figure 
6. We remark that when the cylinder’s surface obeys the SIBC, gi, 2 (D) 
simplify to 


1 e-W* 

9i(D) - 9(D) = -j= J_^ ^(r) _ q Wt(T) dT ^ 

g 2 (D) -* mQg(D) (54) 

In the case of slowly attenuating creeping and/or surface waves (42) 
is not adequate and an improved result may be obtained by adding (42) 
and (43) with t > irb (in case of a circular cylinder of radius b). Clearly, 
the addition of (43) corresponds to the contribution of the creeping wave 
that has travelled the minimum distance on the cylinder’s surface to reach 
Pn- The contribution of those creeping waves that travelled more than 
once around the cylinder is given by u 2 and could be added to if greater 
accuracy is required. 

The functions g\(D), g 2 (D) and G(x,q) are Fock-type integrals that 
are formally the Fourier transform of a slowly varying factor consisting of a 
quotient of terms containing Airy functions and their derivatives. A compu- 
tationally efficient scheme was proposed by Pearson [16] for the evaluation 
of these integrals. The scheme is an extention of the Fourier trapezoidal 
rule devised by Tuck [17] to treat the rotated-ray exponential behavior oc- 
curing in the integrals. The sampling frequency used in the computation is 
dictated by the slowly varying Airy-function factor in the integrands. Suf- 
ficiently accurate results have been obtained for both lossy and non-lossy 
coatings on a perfectly conducting cylinder using this scheme. 
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VII Numerical Results 


The UTD expressions derived in the previous sections provide a complete 
set of equations for the computation of the total field in all regions of 
interest. In this section, we present some calculated data which validate 
the accuracy of the derived expressions by comparison with data based on 
the moment method and eigenfunction solutions. 

In figure 8 the eigenfunction solutions based on the GIBC and SIBC 
simulations are compared with the exact for a coated cylinder and this 
clearly demonstrates the improved simulation (with respect to the standard 
impedance boundary conditions - SIBC) achieved with the second order 
GIBC. To show the validity of the UTD solution in the case of a convex 
cylinder, a special case of an elliptical cylinder (see fig. 7) is considered 
in figure 9. Data based on the moment method are compared with those 
obtained from the UTD solution in conjunction with the second order low 
and high contrast boundary conditions. 

Figure 10 verifies the asymptotic solution developed for the field point 
in the close vicinity of a convex cylinder. We remark, however, that the 
approximations used for the Hankel functions in the derivation of (42) and 
(43) become less accurate for some values of e r and fi T associated with 
lossless coatings, and this can be avoided by using more accurate approxi- 
mations for the Hankel functions. Finally, figure 11 demonstrates the use of 
GIBC in simulating multilayer coatings by simply redefining the material 
constants a m and a' m as discussed in [10, 13]. 

A difficulty in implementing the expressions derived in this paper was 
the evaluation of the Fock-type integrals G(x,g), gi(D) and g 2 (D) as well 
as the determination of the zeros corresponding to (21). The Fock-type 
integrals were evaluated by employing the method described in [16] and 
the zeros of (21) were determined using the routine given in [20]. 

VIII Summary 

Rigorous ray solutions of the scattered fields were presented for a coated 
convex cylinder. These were developed in the context of the uniform ge- 
ometrical theory of diffraction and specific expressions were given for the 
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scattered fields in the lit, shadow and transition regions as well as for obser- 
vations in the near vicinity of the cylinder. That is, UTD expressions were 
derived for all regions exterior to the coated cylinder. These are suited for 
engineering computations and axe given in terms of the generalized Pekeris 
or Fock- type functions whose evaluation was efficiently performed via the 
Fourier Trapezoidal rule suggested by Pearson [16]. 

In comparison to the solution given by Kim and Wang [8], the ray rep- 
resentations given here are based on a second order generalized impedance 
boundary condition which permits the simulation of thin multilayered coat- 
ing as demonstrated in the included examples. Also, in our implementation 
of the transition fields we employed a rigorous rather than a heuristic eval- 
uation of the Fock- type integrals. Further, we have presented accurate field 
representations for observations on or near the vicinity of the coated cylin- 
der and these can also be used for computing the radiated fields by a source 
or an aperture on the surface of the convex cylinder. 
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Appendix 

Consider an elliptical cylinder as an example of a convex geometry shown in 
figure 7. The geometrical and field parameters associated with this surface 
axe listed below: 

Incident field 

At the observation point (p, <f>) a unit amplitude plane wave incident 
at an angle <f>, is specified by (phase of all fields will be referenced to the 
origin) 


ti'(p, <f>) = “•(*-*•) 


(A.l) 
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Geometry 

The surface of an elliptical cylinder is described by a parameterized 
vector as 


r(j) = a cos 7 x + b sin 7 y (A.2) 

where the parameter 7(0 < 7 < 2 tt) is related to the angle <f> via the relation 


tan <f> = - tan 7 

a 


(A.Z) 


Also, a and b are the half lengths of the major and minor axis of the 
elliptical cylinder, respectively. The radius of curvature at any point 7 on 
the elliptical surface is 

(a 2 sin 2 7 + b 7 cos 2 7) 3 ^ 2 


Pg( 7 ) = 


ab 


and the unit tangent and normal to the surface are given by 

r'(7) — a sin jx + b cos 7 y 


t = 


l r '( 7 )l yj a 2 sin 2 7 + ft 2 cos 2 7 

b cos 7 x + a sin 7 y 


n = t x z = 


(AA) 

(A. 5 ) 

(A.6) 


\J a 2 sin 2 7 + 6 2 cos 2 7 
Location of the reflection point Qr 

The direction of the incident ray at the point of reflection Qr is 

s' = — cos 4>iX — sin^ij/ ( A.l ) 

Let Qr be the reflection point corresponding to 7 = 7 r. The direction of 
the reflected ray is then given by 

(p cos <f> — a cos 7 r)x + (p sin ^ — isin 7 R)y 


s = 


(A.8) 


yj{p cos <f> — a cos 7/j) 2 + (p sin <f> — b sin 771 ) 2 

Enforcement of the reflection law n ■ s' = -h-s r at Qr then yields the 
transcendental equation 

p(6 cos <f> cos 7/j + a sin <j> sin 7 r) — ab 


\fp 2 + a 2 cos 2 7 r + b 2 sin 7 r — 2 p(a cos <f> cos 7/1 + 6 sin <f> sin 7/?) 


(A. 9 ) 


= b cos <f>i cos 7/j + a sin 4 >i sin 7/* 
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to be solved for 7 r for each observation angle <j>. The solution must of 
course be done numerically using an iterative routine such as the Newton- 
Raphson’s method. Given 7 r the incident field at Qr 


u‘(Qr) = e‘ U " 


where 


with 


£r = yj a 2 cos 2 ~r + b 2 sin 2 7 r cos (<f> r - fc) 


. , f b 

<p r = tan < — tan 7 r 


} 


The reflected ray caustic distance at Qr is given by 


(4.10) 


(4.11) 


(4.12) 


where 


- Pg(Qri) cos 0' 
" = 2 


(a 2 sin 2 7/* + b 2 cos 2 Jr ) 3 ^ 2 
PgiQn) = : 7 


(A13) 


(A14) 


cos0‘ = n(Q/i) • s r = 


b cos <j>i cos 7« + a sin <f>i sin 7 r 


yja 1 sin 2 7 ~r + b 2 cos 2 7 r 
with £ in (17), (27) and (28) found from 

£ = \J (p cos 4 > — a cos 7/1) 2 + (p sin 4 > — fcsin 7 a) 2 


(A- 15) 


(A16) 


Attachment and detachment points for the surface diffracted ray 
The attachment point Q\ on the surface of the cylinder shown in 
figure 7 can be determined by enforcing the generalized Fermat principle 
implying that at Q x (a cos 71, b sin 71) 


n(Qi)-r = 0 


(A. 17) 
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this yields 


7 i 



(A 18 ) 


Similarly, the surface detachment point Q^(a cos 72, &sin 72) is determined 
by enforcing 


n((J 2 ) ' 5—0 


(A. 19 ) 


where 


s = 


__ ( p cos <f> — a cos 7 2 )x + ( p sin <f> — 6 sin 7 2 )y 


\J{f> cos 4 > - a cos 7 2 ) 2 + (p sin ^ — 6 sin 7 2 ) 2 

This leads to the transcendental equation 

p(a sin <j) sin 7 2 + b cos <f> cos 72) — ab = 0 

which can be solved numerically for 7 2 . 

Incident field at Qi is 


where 


= e jkTl 

ri = \J a 2 cos 2 71 + b 2 sin 2 71 

/ 1 (b ) 

4 > 1 = tan < — tan7i > 


(A.20) 


(> 1 . 21 ) 


(A.22) 


(A. 23 ) 


(A. 24 ) 


Also, the parameters x and t in ( 38 ) and ( 39 ) are, respectively, given by 

ffca 2 6 2 ) 1/3 /"» dj 


x = \— r -4 / (A. 25 ) 

l ^ J * / 'n y a 2 sin 2 7 + 6 2 cos 2 7 

< = F ™ \J a 2 sin 2 7 + ft 2 cos 2 7^7 (A. 26 ) 

The integrals in (A. 25 ) and (A. 26 ) must be computed numerically. 
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a smooth convex cylinder. 



Fig. 4. Reflected ray path. 



Fig. 5. Propagation ray paths in the shadow region. 



Fig. 6. Projection of the field point P in the direction normal to the 
surface at Pyv- 





Fig. 7. Elliptical cylinder configuration. 
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Fig. 8. Bistatic H-polarization scattering pattern for a circular cylinder 
having b = 3A, (a) e r = 4, p r = 1, 6 = 0.07A, p — 5A (b) e, = 8, 
Pr = 1, 6 = 0.2A, p = 3.05A. 
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Fig. 9. Bistatic scattering pattern for an elliptical cylinder having a = 2A, 
b = 1A, p = 5A, <f>i = 0 (a) E-polarization, t T = 4, /i r = 1 , 
6 = 0.07A (b) H-polarization, e r = 8, /i r = 1, 6 = 0.2A. 
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Fig. 1 1 . Bistatic scattering pattern of a three-layer coated circular cylinder 
having b = 3A, & = 0, = 3->0.1, e* = 4-j0.3, ^ = 7-; 1.5, 

Mn = Mr2 - Mr3 = 1, <*i = 0.01A, 6 * = 0.02A, 63 = 0.03A (a) 
H-polarization, p = 5A (b) E- polarization, p = 3.05A. 




